Magnetic field and contact resistance dependence of 
non-local charge imbalance 

A Kleine 1 , A Baumgartner 1 , J Trbovic 1 , D S Golubev 2 , A D 
Zaikin 2 ' 3 and C Schonenberger 1 

1 Institute of Physics, University of Basel, Klingclbcrgstrasse 82, 4056 Basel, 
Switzerland 

2 Institute for Nanotechnology, Karlsruhe Institute of Technology (KIT), 76021 
Karlsruhe, Germany 

3 I E Tamm Department of Theoretical Physics, P N Lebedev Physics Institute, 
119991 Moscow, Russia 

E-mail: andreas . baumgartner@unibas . ch 

Abstract. Crossed Andreev reflection (CAR) in metallic nanostructures, a possible 
basis for solid-state electron entangler devices, is usually investigated by detecting 
non-local voltages in multi-terminal superconductor/normal metal devices. This 
task is difficult because other subgap processes may mask the effects of CAR. One 
of these processes is the generation of charge imbalance (CI) and the diffusion 
of non-equilibrium quasi-particles in the superconductor. Here we demonstrate a 
characteristic dependence of non-local CI on a magnetic field applied parallel to the 
superconducting wire, which can be understood by a generalization of the standard 
description of CI to non-local experiments. These results can be used to distinguish 
CAR and CI and to extract CI relaxation times in superconducting nanostructures. 
In addition, we investigate the dependence of non-local CI on the resistance of the 
injector and detector contacts and demonstrate a quantitative agreement with a recent 
theory using only material and junction characteristics extracted from separate direct 
measurements. 
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1. Introduction 

Transport phenomena in superconductors have been heavily investigated for many years. 
Recently, sub-micron scaled multi-terminal normal metal/superconductor (NS) devices 
came to renewed attention of both, experimental [U, El El HI [5j El [7] and theoretical 
0, El QUI EEH E21 [131 E] investigations, mainly because of a quantum mechanically non- 
local process called crossed Andreev reflection (CAR). In local Andreev reflection, an 
electron impinging on the NS interface from the N side can enter the superconductor 
at subgap energies only by forming a Cooper pair (CP) with a second electron. To 
maintain momentum conservation a hole is retro-reflected into the N contact. As 
depicted in figure la, in CAR the hole is reflected into a second contact N2 separated by 
less than the superconducting coherence length from the first contact Nl. The inverse 
process can be described as splitting of a Cooper pair with the electrons entering two 
separate contacts. Since ideally the two electrons retain their correlations from the 
superconductor, this process might provide a source of spin-entangled electron pairs in 
a solid-state environment [BJ [91 [11] . 

Cooper pair splitting has been demonstrated very recently in a double quantum 
dot system [To] , where the strong Coulomb interaction in the quantum dots leads to a 
suppression of other competing transport mechanisms. In contrast, the experimental 
observation of CAR in metallic structures is more challenging because of other subgap 
processes of similar probability, namely elastic co-tunnelling (EC) and non-local charge 
imbalance (CI). EC is the direct electron transfer from one contact to the other via 
virtual intermediate states in the superconductor. It has been shown theoretically that 
without interactions and to the lowest order in the tunnelling rate, EC exactly cancels 
the effects of CAR. Several recent publications [2j [3j EJ EJ [7] show experimentally that 
this cancellation can be lifted. Various models have been proposed to understand these 
findings [T2 l [TBI [Ti] . for example dynamical Coulomb blockade (DCB) and interactions 
via electromagnetic excitations of the environment [13] . Here we present data that 
support our previous interpretation in terms of DCB [7j. In the main part of the paper, 
however, we focus on charge imbalance (CI), i.e. on non-equilibrium effects in the 
superconductor. It is well-established that at finite bias charged quasi-particles (QPs) 
are excited above the gap in the superconductor by tunnelling injection of electrons, 
which leads to CI in the superconductor [16], see figure lb. These QPs diffuse away 
from the injector contact and lead to a potential gradient in the superconductor, see 
figure lc. Distinguishing CAR from CI and EC is essential for the investigation of CAR 
and can be achieved by ferromagnetic contacts [2], the sign of the non-local signal [Hl[7], 
or the characteristic temperature [7J or distance dependence [21 [3], E]. Usually, several of 
these criteria have to be employed to discern the subgap processes. We show that non- 
local CI exhibits a strong characteristic dependence on a magnetic field applied parallel 
to the axis of the superconducting wire. This dependence is due to the reduction of the 
pair breaking time, well-described by the conventional CI description [16], [19] adapted to 
our non-local measurements. In addition, we study the contact resistance dependence of 
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non-local CI and compare the results to a recent theoretical model of subgap processes 
in superconductors [H]. 



2. Sample characterization and non-local resistance 

The sample was fabricated on a thermally oxidized silicon substrate by electron beam 
lithography and UHV angle evaporation. It consists of a 250 nm wide aluminium wire 
with several palladium contacts fabricated on top of the tunnel barrier prepared by 
oxidation of the aluminum for 15 minutes in a 0.1 mbar oxygen atmosphere. The 
structure is shown in figure Id. The aluminum is about 45 nm thick and has a resistivity 
of 7.0 fiQcm, a mean-free path of 8.5 nm and a 'dirty limit' coherence length of 110 nm. 
In direct transport experiments we find a critical temperature of 1.3 K and a critical 
magnetic field of 240 mT (at T = 0.24 K and with the field applied parallel to the 
superconducting wire). We concentrate on data obtained using the three contacts N4, 
N5 and N6, with the conductance vs. bias characteristics shown in figure le. While the 
junctions N5 and N6 have similar conductances, N4 exhibits roughly five times larger 
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Figure 1. (Color online) Schematic of the non-local processes a) crossed Andreev 
reflection and b) charge imbalance, c) schematic of the spatial distribution of quasi- 
particlc density Q*(x) in a non-equilibrium superconductor device, d) SEM image of 
the sample, e) Barrier characteristics of barriers N4, N5 and N6. The gray dashed 
lines are BTK fits with an additional broadening parameter, f) Normal state junction 
conductance of N6 at base temperature at a magnetic field of 0.5 T> H c , including 
the finite resistance of the Al leads. The lines without symbols show the theoretical 
zero-bias anomaly due to dynamical Coulomb blockade at the indicated temperatures. 
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Table 1. Contact properties. R: barrier resistance at T = 1.5 K; A: contact area; 
RA: resistance-area product; Z: BTK barrier strength; (?ns = Gns/Gnn: normalized 
tunnelling conductance at zero bias and T = 1.1 K. 
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values. The characteristics can be fitted to the BTK model [T7_], with an additional 
broadening parameter and the barrier strength Z used as fit parameters (for N4 we also 
slightly adjusted the 1.5 K normal state resistance R). The superconducting energy gap 
A = 190 yueV and the temperature T = 0.24 K are not used for fitting. Z and R are 
given in Table 1, together with the resistance-area product RA of each junction. All 
barriers are in the regime of moderate transmission, ideal for the observation of CAR 
in our structures [TJ. Figure If shows the normal state resistance of contact N6 at base 
temperature and a magnetic field H = 0.5 T> H c applied parallel to the superconductor 
wire and will be discussed below. The edge-to-edge distances between the contacts are 
e#45 = 105 nm, d^ = 140 nm and d^ = 375 nm. The measurements were performed in a 
3 He cryostat with a base temperature of 0.24 K. 

We measure the non-local resistance as shown schematically in figures lc and Id: 
an ac modulated voltage Ud c + £4c with L7 ac = 12/iVrms is applied to a normal metal 
contact (injector), e.g. N4 in figure Id, which leads to a local current I^c + lac through 
the superconducting wire to ground via SI. We record the ac- modulation U n \ of the non- 
local voltage between the normal contact N5 (detector, not in the current path) and 
the superconductor at S2, more than 10//m away from the normal contacts, by lock-in 
techniques at a modulation frequency of ~ 10 Hz. The non-local differential resistance 
can be approximated by R n \ = U n \j J ac . We expect R n \ < for CAR and R n \ > for EC 
and CI (TJ. 

On the present device we obtain similar R n \ vs. traces as reported previously 
[TJ. As an example, a series of curves at different temperatures is shown for the injector- 
detector pair (N6, N5) in figure 2a. At base temperature R n \ exhibits a local maximum 
at zero bias and becomes negative at finite subgap voltages. R n \ is positive for bias 
potentials eU > A. These characteristics remain unchanged up to T = 0.5 K. At higher 
temperatures the subgap feature is washed out and the signal becomes positive for all 
voltages. Similar curves are obtained for the contact pair (N4, N5), except that no 
zero-bias maximum develops (not shown). We attribute a negative non-local resistance 
to CAR and a positive sign to EC or CI. The development of the local maximum at low 
temperatures for the data with the less transparent injector N6 is probably due to a 
shift in the relative strength of CAR and EC, caused by enhanced dynamical Coulomb 
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To further support this interpretation we report here a zero-bias anomaly in the 
conductance of the junctions N5 and N6, consistent with dynamical Coulomb blockade. 
As an example, we plot the conductance of junction N6 at base temperature in figure If. 
Here, superconductivity is suppressed by an external magnetic field of 0.5 T, to avoid 
features due to the superconductor gap. We find a weak, but clear reduction of the 
junction conductance at zero bias. We also plot the theoretically expected response due 
to dynamical Coulomb blockade for a junction capacitance of 0.5 fF (estimated from a 
plate capacitor model) with a resistive environment [18]. The only fit parameter is the 
environmental resistance, which we adjusted to 105 Q. This value might be explained by 
the transmission-line response of the cryostat cables. The zero-bias anomaly is strongly 
suppressed at higher temperatures. We observe a similar feature in the conductance of 
N5, but a constant conductance for N4. When the latter is used as injector we observe 
no maximum at zero bias in the non-local signal. 

3. Magnetic field dependence of CI 

Now we discuss the magnetic field dependence of the subgap signals. From direct 
conductance measurements of the aluminum strip we infer a critical field of H c m 
240 mT. In figure 2b, a series of R n \ vs. U& c curves of the contact pair (N6, N5) is 
shown for different magnetic fields applied parallel to the Al wire. At bias potentials 
larger than the superconductor energy gap the non-local resistance is strongly reduced 
already for relatively small fields, H « H c . In contrast to this, the subgap features 
are not or only very weakly affected. At H > H c , e.g. for H = 500 mT, the non-local 
response becomes zero for all voltages. 

These findings can be seen clearer in magnetic field traces for the contact pair (N4, 
N5), which exhibits an appreciable non-local resistance at base temperature and zero 
bias, see figure 2c. At T = 0.24 K the signal is negative and we tentatively attribute 
it to CAR. We find that this signal is independent of H for H « H c . Only near H c 
the signal increases and reaches positive values. We note that the signal does not reach 
zero exactly for H > H c , possibly due to a small cross talk current in the device, though 
the signal does not change upon doubling of the measurement frequency. In contrast 
to measurements at base temperature, R n i is positive for T = 1.1 K and H = and is 
reduced strongly for small fields. Near H c the signal increases strongly. For H > H c the 
signal is essentially zero. We obtain a similar magnetic field trace when we apply a finite 
bias to the injector, as shown in figure 2d for U^c = 0.2 mV and U^ c = 0.3 mV at base 
temperature, with critical fields of 205 mT and 190 mT, respectively, and an enhanced 
signal strength due to the larger injected QP currents. 

Our magnetic field data on CAR are in stark contrast to previous experiments [3], 
where CAR was found to be strongly suppressed by a magnetic field parallel to the 
superconducting film. This inconsistent behaviour might be related to the different 
coherence lengths in the materials used for the superconductor wire and normal metal 
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Figure 2. (Color online) Non-local resistance vs. bias for a) a series of temperatures 
and b) a series of magnetic fields parallel to the Al strip for the contact pair (N6, 
N5). c) Magnetic field traces of R n \ at zero bias at T = 0.24K and T = 1.1 K for the 
contact pair (N4, N5). The green line is a fit as discussed in the text, d) Magnetic 
field dependence of R n \ at finite bias and base temperature for the contact pair (N4, 
N5). 



contacts, or the different sample geometries. However, we note that in contrast to our 
previous experiments [7J, the out-of-phase part of the low-temperature signal in figure 2c 
has a similar amplitude as the in-phase signal, so that our results on CAR have to be 
considered as preliminary. 

In contrast to CAR, the magnetic field dependence of CI can be well understood by 
adapting the standard CI description [TB] to non-local experiments: we assume that the 
injected quasi-particles (QPs) diffuse isotropically in the Al wire (ideally, the electric 
field in the superconductor is zero) and that the QP concentration near the detector 
contact determines the induced voltage in a similar way as in local measurements. This 
leads to the following expression for the non-local resistance: 

r , = F *( U dc)r Q * - N /a q * m 

Unl 2e 2 M^NS,det(0) ' 1 J 

where tq* denotes the CI relaxation time, F*(Ua c ) the fraction of current that is carried 
by the injected QPs at a given bias, N the density of states of the normal state Al, 
Qinj the injection volume, for which we insert the superconductor volume below the 
injector contact. gNS,det(0) is the zero bias detector conductance at the temperature 
of the experiment, normalized to the conductance of the junction in the normal state, 
and x is the distance between injector and detector. The charge imbalance relaxation 
length Aq* is the distance over which the QP distribution function in a superconductor 



relaxes to its equilibrium value. It is related to tq* by Aq* + tJDtq*, with D the electron 
diffusion constant in Al. Several processes influence tq*, like inelastic electron-phonon 
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scattering, pair breaking by spin-active impurities, orbital motion, a finite supercurrent 
or inhomogeneities in the order parameter A. Generally, the QP life time reads [T9l [20] 



Around zero bias the injector current is small and we can ignore effects due to 
pair breaking of the supercurrent. For a homogeneous gap parameter this leads to 
T = TcT 1 + (2te) -1 with rg the orbital pair breaking and te the inelastic scattering 
time. For the temperature dependence of the energy gap we use the standard 



BCS expression A(T) = 1.74A(0Wl — T/T c and for the magnetic field dependence 



A(T, H) = A(T)yJl -H 2 /H C (T) 2 with H C (T) = H C (0)(1 - T 2 /T 2 ). The orbital pair 
breaking time as a function of H is given by rs = ft/A(0, 0) • H 2 /H 2 |20j . 

At zero bias the only fit parameters are F* and te and we obtain F* = 0.05 and 
te = 0.25 ns, corresponding to A* = 1.0 fim, for the fit shown in figure 2c. As expected, 
F* << 1 since the total current at zero bias is dominated by Andreev reflection. The 
relaxation time is considerably smaller than reported for thick (~ 12 ns) and thin films 
(~ 4 ns) [20], but consistent with the reduced charge imbalance length found in other 
experiments on superconducting Al wires [H [5]. The reduction of te for thin films was 
attributed to enhanced electron-electron scattering for films of thickness smaller than 



yhD/k B T [20], which we estimate to ~ 160 nm at T = 1.1 K for our structures. We 
expect a similar additional suppression in thin wires. 

The finite bias experiments shown in figure 2d were performed at base temperature, 
which leads to a larger energy gap and a characteristic increase of CI at a larger field 
compared to figure 2c. In addition, the finite supercurrent at an increased bias leads to 
a reduction of the energy gap and we expect a reduction of the QP relaxation time due 
to additional pair breaking. The latter effect is strong only for large currents near H c , 
so that we can attribute the increase in R n \ with bias at H = to a change of F* : for a 
larger bias a larger fraction of the total current is generated by QP injection. We note 
that R n \ turns negative when the injected current is comparable to the critical current 
for given external parameters. We observe this effect in figure 2d for H > H c , as well 
as in the bias sweeps at elevated temperatures shown below in figure 3a. We attribute 
this phenomenon to the initial stages of the breakdown of superconductivity, possibly 
at phase-slip centers, which renders the current inhomogeneous. 

4. Contact resistance dependence of charge imbalance 

The previous section shows that the QP diffusion is well described by the standard 
theory. We now turn our attention to the contact resistance dependence of the CI 
induced non-local resistance. This is an important question since CAR and EC are 
strongly affected by the junction characteristics [7]. 

To obtain a reasonable signal strength, we performed the experiments at T = 1.1 K. 
In figure 3a the non-local resistance R n \ is plotted vs. bias for a series of injector-detector 




(2) 
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pairs. The curves with the same injector contact are almost identical and can be clearly 
distinguished from the other curves. The distinction is most obvious if we use N4 as 
injector, with a considerably larger conductance than the other contacts. In contrast, 
using N4 as detector does not lead to a significant change in the non-local resistance, 
as can be seen for example for the pairs (N5, N4) and (N5, N6). A comparison of the 
non-local resistance at zero bias between injector and detector pairs as a function of 
temperature is shown in figure 3b. We note the characteristic shape of the curve that 
allows one to identify CI [H, [5]. We find that near zero bias and for T > 0.5 K the non- 
local resistance becomes the same for all contact pairs and thus seemingly independent of 
the contact resistances. For bias potentials larger than the energy gap the corresponding 
curves are not identical and the signal amplitude depends stronger on the injector than 
on the detector resistance, see figure 3a. 

First we compare our data with the standard theory described by equation 1 and 
investigate the relative change between two curves. For the measurements with the 
same injector, F* and f2 in j are identical and we expect that only the distance between 
the contacts and the normalized conductance ^NS,det are relevant. As an example, using 
the experimentally determined <?NS,det at T = 1.1 K and the CI length obtained in the 
previous section, we expect a factor of ~ 1.3 between data points from the pairs (N5, 
N4) and (N5, N6) at the same bias. In figure 3c these two curves are replotted on a 
smaller scale. We obtain a factor 1.0 ± 0.3 between these two curves at zero bias. The 
standard deviation is even larger for the other pairs. We therefore conclude that the 




Figure 3. (Color online) Non-local resistance a) vs. bias at T = 1.1 K and b) vs. 
temperature at zero bias for the indicated injector-detector pairs, c) and d) show the 
corresponding curves calculated in the model of reference [14j with the experimental 
traces discussed in the text. 
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deviation from the standard description of CI is not significant due to the low signal-to- 
noise ratio. At higher bias the deviation from equation 1 is more pronounced, as shown 
in figure 3a, where the expected non-local resistance at U dc = 0.8 mV is indicated by 
a green and a blue arrow for the pairs (N5, N6) and (N6, N5), based on the values of 
the pairs (N5, N4) and (N6, N4), respectively. However, in this regime an appreciable 
supercurrent is flowing in the device and we do not expect equation 1 to hold exactly. 

Swapping the injector and detector contacts allows one to compare signals 
independent of the contact distance. In this scheme, however, it is necessary to estimate 
F* and the injection volume, each introducing considerable errors. For example, by 
assuming that F* is similar for two contacts (e.g. for large bias F* « 1) one finds a 
factor of ~ 3.5 between the curves with swapped contacts (N4, N6) and (N6, N4). This 
is not supported by our data, see figure 3b, which shows that F* varies between the 
contacts even at zero bias. 

Our experimental results can be understood quantitatively by the model put 
forward in reference [Hj, which we adapt slightly so that the experimental junction 
parameters in table 1 can be used. The non-local resistance of our device in terms of 
conductances reads 

Rn\(U dc ) = G n\{Udc) ^ 

Gdet(0)Gmj(^dc) — Gn^t/dc)' 



where 



GJU a ) = — [ dE 9ai - E ) P rr (4) 

R a J 4fc B Tcosh 2 E - eU « 



2k B T 

(a =inj, det) is the local conductance of the injector and detector and g a (E) are the 
energy dependent spectral conductances of the NS barriers |17j : 

2e(A-\E\)(l + Zl)A* 
9a{ ' £ 2 + (A 2 -£ 2 )(l + 2Z 2 ) 2 
29{\E\ - A){1 + Zj)\E\ 

|£| + V£ 2 - A 2 (1 + 2Z 2 )' { ' 

with 9 the Heaviside function. The general expression for the non-local conductance 
G-a\_(U dc ) is presented in reference [H]. Provided the charge imbalance length Aq» is 
shorter than the length of a superconducting wire (which is the case in our experiments), 
we find 

1 



G n \{U dc 
X / dE 



4e 2 NDSR in] R det 

g inj (E)g det (E) A 2 - E 2 e~ k ^ 
4fc B Tcosh 2 4^ A 2 k(E) 

\E\<A 2k B T y ' 

Ao*e~^ A Q* 



4e 2 NDSR inj R det 



g mi (E)g dct (E) ff 2 -A 2 
4A; B Tcosh 2 ^^ E* ' W 
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where k(E) = y2V A 2 — E 2 / D + 1/Aq* and 5 is the superconducting wire cross- 
section. Note that at zero bias and in the limit A <C k^T, i.e. close to the critical 
temperature or the critical magnetic field, equation [3] reduces to equation [T] with 
E* = f2 in j/2 1 S , A<2». One can see that for G a 3> G n \ the resistances of the NS barriers 
cancel out and do not enter the expression for the non-local resistance, which depends 
on the barrier properties only through the BTK parameters Z a . Furthermore, in both 
limits of high and low temperatures also these factors drop out and one arrives at the 
following simple expressions for the non-local resistance 



where r\ Q:t and rg are the normal state resistances of the superconducting wire segments 
of lengths Aq* and £, respectively. Thus, we find that the non-local resistance in these 
limits only depends on the distance between the junctions and the properties of the 
superconducting wire, but not on the properties of the NS barriers. 

In figure 3c, numerically calculated _R n i curves based on the above equations are 
shown for T — 1.1 K, with material and contact characteristics from the experiments, 
including Aq* = 1 /im from the magnetic field dependence. We stress that no fit 
parameters are used for these plots. We find that the calculations quantitatively 
reproduce our data. The two curves with larger contact separation exhibit a reduced 
non-local resistance compared to the other curves, which could not be resolved in the 
experiment. Figure 3d shows calculated temperature sweeps, which compare very well to 
the experimental curves. For this comparison we have subtracted a small constant from 
the data, since the model does not account for negative non-local resistance, because no 
electron interactions are considered in this model. In this model, the swapping of the 
injector and detector results in identical curves, also in agreement with our data. 

5. Conclusions 

In summary, we report measurements of the non-local resistance in a multi-terminal 
superconductor/normal metal device. We show a zero-bias anomaly at base temperature 
in the normal state of the superconductor consistent with dynamical Coulomb blockade, 
which might be partly responsible for lifting the balance between crossed Andreev 
reflection (CAR) and elastic co-tunnelling. We then focus on charge imbalance (CI) 
and demonstrate a characteristic magnetic field dependence, which provides a novel 
tool to distinguish CAR and CI in future experiments. By comparing our data to 
a generalization of the standard description of CI, we can extract the inelastic quasi 
particle relaxation time and find a reduction relative to literature values probably due 
to the reduced dimensions of the superconductor wire. 

In addition we present a set of experiments where we analyze systematically pairs 
of injector and detector contacts. To generate CI, we chose a rather high temperature, 
at which the superconducting and the normal junction resistances differ only weakly, 
which results in only small differences in the curves. The results suggest that the 




T c - T < T c , 
k B T < A, 



(7) 
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contact resistances influence the non-local CI only weakly. This finding is reproduced in 
a theoretical model adapted to incorporate experimentally accessible parameters. This 
model does not require additional adjustable parameters and agrees quantitatively with 
our data. 
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